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Abstract — In this paper, we investigate a unified linear 
transceiver design with mean-square-error (MSE) as the objective 
function for a wide range of wireless systems. The unified design 
is based on an elegant mathematical programming technology 
namely quadratic matrix programming (QMP). It is revealed 
that for different wireless systems such as multi-cell coordi- 
nation systems, multi-user MIMO systems, MIMO cognitive 
radio systems, amplify-and-forward MIMO relaying systems, the 
MSE minimization beamforming design problems can always be 
solved by solving a number of QMP problems. A comprehensive 
framework on how to solve QMP problems is also given. 

Index Terms — Quadratic matrix programming, transceiver 
design, optimization, MSE. 

I. Introduction 

In order to satisfy the forecasted data rate boost and 
to enable high quality multi-media wireless services, more 
and more available wireless resources are introduced into 
wireless systems, such as temporal, frequency and spatial 
resources. The multi-dimensional wireless resources bring new 
challenges for wireless system designs. To order realize the 
promised performance of these resources, the corresponding 
new technologies will be adopted, such as multiple-carrier 
technology, multiple-antenna technology and so on. Q 

For spatial resource, multiple-input multiple-output 
(MIMO) technology is a great success in both theoretical 
research and practical applications Along with the 

evolvement of wireless systems, MIMO becomes to be 
a fundamental and important ingredient of complicated 
wireless systems e.g., cooperative communication, cognitive 
communication, physical layer security communication, 
network coding based communication and so on. With 
channel state information (CSI), transceiver designs can 
significantly improve the performance of MIMO systems 0. 

Because of a variety of wireless service requirements and 
wireless environments, there are various wireless systems 
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which have totally different network architectures with even 
different wireless interfaces. For all of these wireless systems, 
transceiver design can significantly enhance the whole system 
performance J2J. However, for different wireless systems the 
transceiver design problems admit different signal models, 
different power constraints and different numbers of variables. 
It means that these designs should be investigated case by case. 
From theory research perspective, the theorists want to find a 
unified design which can reveal some nature of the transceiver 
designs. To the best our knowledge, the transceiver designs are 
not unified simultaneously from different performance metrics 
and different systems. 

In the existing works, for unified linear transceiver designs, 
the widely used logic is for a given system linear transceiver 
designs with different performance metrics are unified into 
a general caless of optimization problems with the similar 
mathematical formulations 0, Q. It is well-known that 
there are two guidelines, i.e., using magorization theory 
and weighting operation (TJ). When magorization theory is 
used, the transceiver design logic is to formulate different 
performance metrics as different functions of the diagonal 
elements of the data detection MSE matrix. Furthermore, the 
functions are classified into two main subclasses of function 
i.e., Schur-convex or Schur-concave functions. On the other 
hand, using weighting operation, the different performance 
metrics are taken into account by using different weighting 
matrices. 

In contrast to these works, in this paper, we give a unified 
transceiver design which aims at unifying the linear transceiver 
designs for different wireless systems with the same perfor- 
mance metric minimum mean square error (MMSE). It is 
revealed that for the beamforming designs in different wireless 
systems such as multi-cell coordination beamforming design 
problem, multi-user MIMO beamforming design, cognitive 
MIMO beamforming design, amplify-and-forward MIMO re- 
laying beamforming design and even their corresponding 
robust transceiver designs with randomly distributed channel 
estimation errors and so on, the transceiver design problems 



can be optimized by solving a series of quadratic matrix 
programming (QMP) problems. By the way, we also clarify 
that the QMP problems can always be efficiently solved. 

The following notations are used throughout this paper. 
Boldface lowercase letters denote vectors, while boldface 
uppercase letters denote matrices. The notations Z T , Z* and 
Z H denote the transpose, conjugate and conjugate transpose of 
the matrix Z, respectively and Tr(Z) is the trace of the matrix 
Z. The symbol Im denotes an M x M identity matrix, while 
Om.n denotes an M x N all zero matrix. The notation Z 1//2 is 
the Hermitian square root of the positive semi-definite matrix 
Z, such that Z 1 / 2 Z 1 / 2 = Z and Z 1 / 2 is also a Hermitian 
matrix. The symbol E denotes statistical expectation operation. 

II. Motivations 

At the beginning, we would like to discuss why our 
attention is concentrated on linear minimum mean-square- 
error (LMMSE) transceiver designs in this paper. However 
for ceratin performance metrics such as bit error rate (BER) 
the performance of linear transceivers may be not as good 
as that of their nonlinear counterparts, linear transceivers are 
still preferred by practical wireless systems due to their low 
implement complexity. 

On the other hand, MSE is a widely used performance 
metric of the optimization problems for estimation, detection 
and transceiver design. It should be pointed out that MSE 
acting as performance metric suffers from several inherent 
drawbacks as it is not a terminal performance metric e.g., 
capacity and BER. Roughly speaking, MSE can be seen as an 
approximation of the final performance metrics, although they 
have very closed relationships and particularly, in some special 
cases, they are even equivalent with each other. Furthermore, 
MSE is a statistical average, which cannot reflect the detailed 
information for specific samples. 

Maybe more weaknesses can be found, if we discuss in 
more depth, however MSE is still preferred by the wireless 
researchers. It may be not the best, but it is not the worst. In 
general, tractability is the main advantage of MSE. In com- 
plicated wireless systems, for several terminal performance 
metrics, their formulations may be too complicated to prohibit 
them from application. From engineers' point of view, the case 
where there is a solution is much better than that where there 
is no solution. 

In the following, we show an example to illustrate that 
QMP problems are of great importance in LMMSE transceiver 
designs. Note that this example is discussed in detail in 
our previous work Q. Here, it only provides a prologue 
of our work in this paper. First, we want to highlight that 
the algorithm discussed in this paper is not limited to the 
example, which has a much wider application range. For 
example, the algorithm discussed in this paper can also be 
applied to multi-hop networks. Moreover, different from our 
previous work @), in this paper we provide a comprehensive 
framework for LMMSE transceiver design. We are not limited 
to a specific communication systems. We try to reveal the 
nature of LMMSE transceiver design and answer the questions 



why QMP should be chosen and how to solve the optimization 
problem using QMP. 

A. An example: 

In the example, a dual-hop amplify-and-forward (AF) relay 
network is considered, in which there are multiple source 
nodes, relay nodes and destination nodes. Furthermore, dif- 
ferent source can have different numbers of transmit antennas 
and data stream to transmit. It is denoted that the number 
of transmit antennas of the i th source is iVgj. It is also 
assumed that for each source node there may be more than 
one corresponding destination node. There are also multiple 
relay nodes in the network, and the j th relay has Mrj receive 
antennas and Nrj transmit antennas. At the first hop, the 
source nodes transmit data to the relay nodes. The received 
signal, x 3 , at the j th relay node is 

+ ni d . (1) 

where Sjfc is the data vector transmitted by the i th source 
node to the fc th destination with the covariance matrix R Si) , = 
EjsifeS^,}. When the i th source node does not want to transmit 
signal to fc th destination, Sjfc is a all-zero vector. 

At the source, before transmission the signal is mul- 
tiplied a precoder P^- under transmit power constraint 
J2k Ti"(PifeRs it .Pife) < Ps,i, where P S i is the maximum 
transmit power at the i th source node. The matrix H sr> jj is the 
MIMO channel matrix between the i th source node and the j th 
relay node. Symbol is the additive Gaussian noise with the 
covariance matrix R ni . At the j th relay node, the received 
signal Xj is multiplied by a precoder matrix Fj , under a power 
constraint T^Fj-R^F?) < P r j where R Xj . = E{xjxf} and 
P r j is the maximum transmit power. Then the resulting signal 
is transmitted to the destination. The received signal at the k th 
destination, y k , can be written as 

yfe = } j XH r d,jk'Fj-x.j) + n 2 ,, 

= 2j,[H rd ,j fc Fj^ (H^jjPjfcSjfc)] 

+ /] . [Hrd,jtl?j / (H ar jj) (PzmS/ m ))] 
* — * — *l z — 'm^fc 

+ ^.(Hr^feFj-nij) + n 2 , fc . (2) 

where H r djk is the MIMO channel matrix between the j th 
relay and the fc th destination, and n-2.k is the additive Gaussian 
noise vector at the second hop with covariance matrix R n , h . 

The optimization problem of linear minimum mean-square- 
error (LMMSE) transceiver design can be formulated as (4) 

mm £ fc MSE fc = E{||G fc y fc - [ s J k , • • • , s^ fc ] T || 2 } 

s.t. Tr(F J R x .Ff ) < P rd j e £ r 

^ fc Tr(P. tfc R Sifc Pf fc ) < P Vi ie£ s (3) 

where [sj k , ■ ■ ■ ,sjf k ] T is the desired signal to be recovered 
at the fc th destination, and £ r and £ s denote the set of relay 
nodes and the set of source nodes, respectively. 



The optimization problem |3) is a very general problem 
which includes the following scenarios as its special cases. 

• Multi-user MIMO uplink transceiver design J6). 

• Multi-user MIMO downlink transceiver design |6j. 

• Two-way AF MIMO relaying LMMSE transceiver design 
15) . Two-way AF MIMO relay can be taken as a soft combi- 
nation of uplink and downlink beamforming designs. 

• Multi-cell coordination beamforming design. 

As there are too many variables to be optimized and the 
nonconvex nature of the optimization problem makes the 
problem more complicated, it is difficult to find a closed- 
form globally optimal solution. In order to carry out the 
transceiver designs, several suboptimal solutions are usually 
proposed. Iterative algorithm is one of the most widely used 
and important suboptimal solutions. We admit that iterative al- 
gorithms suffer from some well-known weaknesses. First, the 
final solution is greatly affected by the initial value selection. 
Second, the convergence of an iterative algorithm should be 
guaranteed. If not, the iterative algorithm may be meaningless. 
Third, in general even with proved convergence there is no 
guarantee that the final solution is globally optimal. However, 
iterative algorithms still have two important characteristics 
making them preferable. First, it can be applied to a much wide 
area of transceiver designs ranging from traditional a point- 
to-point system to a distributed network. Second, it can act 
as a performance benchmark for other suboptimal solutions. 
Inspired by these, we focus on iterative algorithms in this 
paper. 

B. Quadratic nature of the LMMSE transceiver designs 

The MSE is an integration over the signals and noises. From 
its name, it is direct and obvious that MSE is a quadratic 
formulation. Moreover, in this paper, we concentrate our 
attention to the case where the variables are matrices, as 
in MIMO systems the variables to be optimized are always 
matrices. Inspired by these facts, a kind of function termed 
as quadratic matrix (QM) functions with a complex matrix 
variable X is defined as 

/,(X) = Tr(D ; X H A ; X) + 27e{Tr(B[ i X)} + q (4) 

where A ; = Af e C nxn , B; G C nxr , q e M, D; = 
Dp e C rxr . In addition, denotes the real part. It can 

be seen that a QM function consists of three terms which 
are second-order term, first-order term and zeroth-order term. 
If the following conditions are satisfied, no matter what the 
system is, the MSE with linear transceiver is a QM function 
with respective to each variable, separately. 

(1) . The considered system is a linear system. The linearity is 
mainly referred to the following two properties: 

(a. 1) The received signal at the destination is a linear 
function of the transmit signal when all variables are fixed. 

(a.2) The received signal at the destination is a linear 
function with respective to each variable when the signal is 
fixed and the other variables are fixed. 

(2) . The desired signals are independent of the noises. It means 
that when the signal vector is denoted by s and the equivalent 



noise vector is v, the following equality must hold 

E{sv H } = 0. (5) 

In addition, the constraints in transceiver design for wireless 
systems are usually QM functions, as the constraints are 
usually related with energy, which are definitely quadratic 
terms, e.g., transmit power, interference to primary users, and 
so on. Therefore, most of LMMSE transceiver designs can 
be iteratively optimized by solving a number of optimiza- 
tion problems consisting of QM functions in both objective 
function and constraint functions. This kind of optimization 
problem is named as quadratic matrix programming (QMP) 
problems. In this following, the properties of QMP will be 
discussed in detail as well as how to solve it. 

III. Fundamentals of QMP 

At the beginning of this section, the definition and properties 
of QMP are investigated. Although in 0, a definition of 
quadratic matrix programming is given, in this paper we first 
revise the definition given in J7] in order to accommodate more 
cases. Our definition is more general and has a wider range 
of applications. A standard QMP problem is defined as 

Type 1 QMP: 

min Tr(DoX H AoX) + 2ft{Tr(B^X)} + c 

s.t. Tr(D 1 X H A l X) + 27e{Tr(Bf X)} + a < 0, i G 1 
Tr(D 2 X H A J X) + 2ft{Tr(BfX)} + Cj = 0,j G £ 
X G C" xr (6) 

where I G {0} UlU£, These assumptions are essential to 
guarantee that the objective function and constraint functions 
are real-valued functions, as it is meaningless to minimize 
a complex-valued function. It is obvious that QMP is a 
special case of quadratically constrained quadratic program- 
ming (QCQP), which is a class of very famous and widely 
used optimization problems. Then it is natural that QMP has 
much better properties than QCQP, which can be exploited to 
solve the involved optimization problems. This is exactly the 
motivation of the research on QMP Q, J8). 
General QMP: 

Based on the properties of Kronecker product and the 
following definitions 



ft. 



Di<g)A ; vec(B z ) 
vec H (B ; ) ci 



(7) 



the optimization problem © is equivalent to 
min Tr(JloZ) 

s.t. Tr(f2iZ) < 0, Tr(J2j-Z) = 
Z = [vcc T (X) l] T [vcc H (X) 1]. 



(8) 



If the constraint Rank(Z) = 1 is relaxed (it is a well-known 
semi-definite relaxation (SDR) (SJ), we have the following 



semi-definite programming (SDP) problem iflOl , which can 
be efficiently solved by interior point polynomial algorithms 



In the remaining part of this paper, we will concentrate our 
attention to a special QMP defined as 



min Tr(O Z) 

s.t. Tr(niZ) < 0, Tr(%Z) = 

[ZJjvjv.+ijivjv.+i = 1, z y 0, (9) 

where Z is a Hermitian matrix. When A/ and D/ are both 
positive semi-definite matrices and there are only inequality 
constraints, the optimization problem (|6]l is a convex opti- 
mization problem. For convex optimization problem, we do 
not need the previous SDR to solve the optimal solutions. In 
the following, two methods for the convex case are proposed. 
Convex QMP: 

When A; and D; are both positive semi-definite matrices 
with inequality constraints, the constraints only exists inequal- 
ity constraints. 
SDP Based Algorithm: 

Using the properties of Kronecker product Tr(AB) = 
vec H (A H )vec(.B), the QM function can be reformulated as 

Tr(Df /2 X H A ; XD, 1/2 ) + 2ft{Tr(Bf X)} + q 



Tr(Df /2 X H Af A ; 2 XD; 1/Z ) + 2ft{Tr(Bf X)} + q 
vec H (X)(D; /2 ® A? XDp ® Af )vec(X) 



1/2, 



27e{vec H (B ; )vcc(X)} + c ( < 0, 



(10) 



based on which and together with Schur complement lemma, 
the optimization problem (01 can be reformulated as the SDP 
problem given at the top of the next page. 

Notice that in our work, the variables are complex matrices. 
For some optimization tool boxes, only real variables are 
permitted. In that case, a minor transformation is needed, 
which is 



v 

a 



>- 



l 2N 



V 

a 



y o 



where v is defined as 



v = [Rcal(v) T Imag(v) 



TiT 



(12) 



(13) 



When Aj and D; are positive definite matrices, the op- 
timization problem can be further transformed into a more 
efficient convex optimization problem termed as second order 
conic programming (SOCP) problems. 
SOCP Based Algorithm: 

Notice the QM functions in both the objective function and 
constraints can be reformulated as 

Tr(Df /2 X H A,XDj /2 ) + 27e{Tr(B[ i X)} + c x 



A, 2 XD; 



+ Cl 



T^A^BjD^Bf) 



(14) 



where || • ||p denotes the Frobenius norm. Therefore, the 
optimization problem (O can be reformulated as a standard 
SOCP problem which is given at the top of the next page. 



Type 2 QMP: 

min Tr(X H A X) + 2ft{Tr(B"X)} + co 

s.t. Tr(X H A,X) + 2ft{Tr(Bf X)} +c,<0,iel 
Tr(X H A J X) + 2-£{Tr(B?X)} + c = 0,j e £ 
X G C nxr (16) 

which has much better properties than the general Type-1 
QMP J7) . For notational simplicity, in the remaining part of 
the paper, the T-2-QMP problems are referred to as the Type 
2 QMP problem. 

IV. Properties of T-2-QMP 

A. T-2-QMP without Constraints 

At the first glance, we discuss the case without constraint 
which reads as 



min Tr(X M A X) 



27e{Tr(B^X)} 



Co 



(17) 



where Ao > 0. This case corresponds to linear minimum 
mean square error (LMMSE) equalizer design, which is also 
named as LMMSE estimator design. It is obvious that the 
optimization problem (fTTI i is convex and then the optimal 
solution is exactly the solution satisfying its differentiation 
equals to 0, i.e., AoX = Bo- Specifically, the optimal solution 
has the following closed-form solution 



L opt 



A^Bq. 



(18) 



Moreover, considering weighted MSE minimization, the 
optimization problem becomes to be 

minx Tr(WX H A X) +27e{Tr(W H BHx)} + c (19) 

where W >z is a weighting matrix. Following the same 
logic as previously discussed, we have the optimal solution 
must satisfy 



A XW = B(,W 



(20) 



Because W can be ill-rank, the optimal solution is not unique. 
It is obvious that X opt = A^Bo satisfying the previous 
condition, which is exactly the optimal solution. 
Conclusion 1: Without constraints, the optimal solution X opt 
has a closed form. Notice that X^ pt is just the Wiener filter. 
It is well-known for a linear system with Gaussian noise, 
LMMSE equalizer is exactly the optimal equalizer ifTTl . 

B. T-2-QMP with One Constraint 

In this section, we focus on the case where there is only 
one constraint in a QMP problem. This case corresponds to 
the scenario there is only one transmit power constraint. Here 
we focus on the following T-2-QMP problem 



min Tr(X H A X) + 27e{Tr(B X)} + c 



x 



s.t. Tr(X H AiX) < P, 



(21) 



s.t. 



_ ((D 2 ^A^ve^X)) 1 



((D/ ®Af)vec(X)) 1 



(D| ® A|)vec(X) 
-27e(vec H (B )vcc(X)) +t 

(Df ®A?)vec(X) 
-2ft(vec H (B 8 )vec(X)) - c, 



^ 



h 



(11) 



min t 
PkJ 

s.t. 



A 2 XT) 2 — A 



A„ 2 XD, 2 



B D 



15 I) 



< t 



< 



(15) 



where Ai > 0. Without loss of generality, it is assumed 
that the feasible set of the considered optimization problem 
is not empty. In this scenario, solving an unknown matrix 
variable can be reduced to solving an unknown scalar variable. 
Therefore, the computation dimensionality and complexity are 
significantly reduced. In this following, we will clarify this 
step by step. 

For constrained optimization problems, when regularity 
conditions are satisfied, Karush-Kuhn-Tucker (KKT) are the 
necessary conditions for the optimal solutions and then they 
can provide useful information to derive the optimal solutions. 
For the optimization problem (l2T~b with one constraint, the 
regularity condition termed as linear independence of con- 
straint qualification (LICQ) can be easily proved to hold 
and thus KKT conditions are the necessary conditions. The 
corresponding Lagrange function of the optimization problem 
(f2Tb is expressed as 

£(X) = Tr(X H A X) + 27e{Tr(B X)} + c 

+ /i(Tr(X H A 1 X)-P), (22) 

where (j, > is the Lagrange multiplier. Based on ( l22l . 
the KKT conditions of the optimization problem (f2TT > can be 
derived to be J9) 

(A + /iAi)X = B£ (23) 

A i(Tr(X H A 1 X) - P) = (24) 

Tr(X H A!X) < P (25) 

H > (26) 

In this case, the optimal solution has the following semi- 
closed-form solution 

X=(Ao+/iAi)- 1 B^ 1 (27) 

in which the only unknown variable is a scalar Lagrange 
multiplier. Substituting d27l i into the constraint of (BTT i. we 



have 

Tr(X H A!X) 

=Tr(B (A + ^A 1 )" 1 A 1 (Ao + nA^Bf) 

=Tr(B A^(A^A A^ + A J)- 2 A^B«) 

=<?(/*)• (28) 

In Appendix A, it has been proved that g(fi) is a decreasing 
function with respective to /i, and its value can be computed 
by using a simple one-dimensional search such as bisection 
search. Based on this conclusion and the KKT conditions given 
previously, the value of /i can be computed to be 

,= {° , if3( ° )£P . (29) 

I Solve g(fi) = P Otherwise 

It can be seen that the solution satisfying KKT conditions 
is unique. Notice that the KKT conditions are the necessary 
conditions for the optimal solutions. As a result, the unique 
solution satisfying the KKT conditions is exactly the optimal 
solution. This is of great importance, because the unknown 
variable is simplified from a matrix to a scalar. In other 
words, the number of variables is significantly reduced and the 
corresponding computational complexity is also significantly 
reduced. 

Conclusion 2: With only one constraint, the T-2-QMP prob- 
lem has a semi-closed-form solution with an unknown scalar 
variable. The unknown scalar variable can be efficiently com- 
puted using one dimensional search. 

Remark: Notice that in Boyd's book J9), it never states 
that the KKT conditions are the necessary conditions for the 
optimal solutions without any prior conditions. It is possible 
that KKT conditions are not the necessary conditions. 

C. T-2-QMP with more than one constraint 

For T-2-QMP problems with more than one constraint, 
solving the optimization problems must also rely on numerical 
algorithms such as interior point algorithms. As a T-2-QMP 



problems has much better structures compared to the general 
QMP problem discussed in Section [TTTJ it exhibits stronger 
convexity property which can be exploit to efficiently solve 
the optimization problem. As discussed in (8), the original 
optimization problem is first transformed into its homogenized 
problem which can be efficiently solved. First, the homoge- 
nized QM function of the QM function defined previously is 
denoted by / 2 H 

/; H (Y; Z) =Tr(Y H A;Y) + 2ft{Tr(Z H Bf Y)} 
+ Q / 7 Tr(Z H Z). 

Then introducing the following operators, 



(30) 



M,(/i) = 



A, B, 

' 7 r 



(31) 



The homogenized optimization problem of ( fT6b is formulated 

as 

min Tr(M(/ )[Y;Z][Y;Z] H ) 
s.t. Tr(M(/i)[Y;Z][Y;Z] H ) <a u iel 
Tr(M(/,)[Y;Z][Y;Z] H )=a„. ? G£ 
Z H Z = I r Y 6 C" xr . (32) 

Notice that the optimal solution of ([Tol l X opt equals to 

X op t = Y opt ZQ pt . (33) 

Defining U = [Y;Z][Y;Z] H , after relaxing the rank con- 
straint on U, we have the following optimization problem 

min Tr(M(/ )U) 

s.t. Tr(M(/i)U) < on,i el 
Tr(M(f j )V) = a j ,jeS 



+ l:n+r,n+l:n+r 



i, 



u y o. 



(34) 



To recover X from U, a rank reduction based algorithm has 
been given in Q. When the number of the constraints are less 
than 2r, this relaxation is tight (8). It is concluded that T-2- 
QMP has a much stronger convexity property than the general 
QMP. 

Remark: From the practical viewpoint, due to the limited 
length of training sequences and the time varying nature 
of wireless channels, channel estimation errors are always 
inevitable. When channel errors are taken into account, the 
channel state information can be written as H; = H; + AH; 
where fi/ is the estimated H; and AH/ is the corresponding 
channel estimation error, respectively. The channel estimation 
errors are independent of the signal and the noise. Notice that 
matrix expectations keep the quadratic nature of the original 
QMP problems lfl~3Tl and therefore QMP still works for robust 
beamforming design ||T2l . 

V. Conclusions 

In this paper, we discussed linear transceiver designs with 
MSE as the performance criterion. Different from the previous 
existing works, the transceiver designs were understood from a 



unified optimization problem named QMP problems for vari- 
ous wireless systems. The QMP based transceiver design algo- 
rithm discussed in this paper can be applied to a wide range of 
problems such as multi-cell coordination beamforming design, 
multi-user MIMO beamforming design, cognitive radio MIMO 
beamfoming design, cooperative network beamforming design 
and even their corresponding robust designs with randomly 
distributed channel estimation errors. Meanwhile, the elegant 
properties of QMP problems were also discussed in detail in 
our work. 

Appendix A 

MONOTONICITY OF g(fl) 

In order to prove that <?(/z) is a monotonically decreasing 
function with respect to /.t, we assume that [i\ > fi2 and 
in the following we will prove that g(pti) < p(/i2)When 
fii > based on the fact that A is a positive semi- 



definite matrix we have lfl4l (A x 1 ^ 2 AqA 1 



-1/2 



(a~ 1/2 a a- 



-1/2 



Mil) > 



/i2i), based on which, taking inversion of 



both sides, the following inequality holds (A 1 AqA 1 + 
fiil)- 1 < (A~ 1/2 A A~ 1/2 + Therefore, it can be 

concluded that 

TrlBoA-^A-^AoAr 1 / 2 + Mi^A^X] 
< Tr[B A- 1/2 (A- 1/2 A A- 1/2 H 
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